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Abstract. Infrared spectra of NGC 7027 and NGC 6543 rang- 
ing from 2.4 to 45 /xm were obtained with the Short Wavelength 
Spectrometer on board the Infrared Space Observatory. A first 
analysis of these spectra, with the aid of photo-ionization mod- 
els, is presented. 

We report the first detection of the [Arvi] 4.53 fim and 
[Nevi] 7.65 ^m lines in the spectrum of NGC 7027. When 
compared with older observations it is clear that the [Ar VI] line 
and possibly also other lines have increased in strength since 
1981. We argue that a likely explanation for this variability is 
a change in the spectral energy distribution of the central star, 
possibly an increase in effective temperature. However, this 
result needs to be confirmed by further observations. 

We also report a non-detection of the [Orv] 25.9 /im line 
and the first detection of the [Nam] 7.32 /im line in the spectrum 
of NGC 6543. The non-detection is not expected based on a 
blackbody approximation for the spectrum of the central star 
The ionization threshold for 0^+ is just beyond the He II limit, 
and the absence of this line shows that the stellar flux drops at 
least by a factor 350 at the He II limit. ModeUng the [O IV] Une 
may prove to be a valuable test for atmosphere models. 

Key words: Stars: evolution - Stars: atmospheres - planetary 
nebulae: individual - Infrared: ISM: lines and bands 



1. Introduction 

Both the planetary nebulae NGC 7027 and NGC 6543 are very 
well studied nebulae. However, a complete study of the infrared 
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spectrum of these sources has not been possible until now. In 
view of their scientific interest and their brightness they were 
obvious candidates for early targets to be observed by the Short 
Wavelength Spectrometer (SWS) on board the Infrared Space 
Observatory (ISO). In this paper we will present selected ionic 
emission lines of both SWS spectra. We will also present pre- 
liminary photo-ionization models and we will discuss possible 
interpretations of the reported features. 

2. The ISO SWS observations 

The SWS spectra of NGC 7027 and NGC 6543 were obtained 
during the Performance Verification phase, in the 24th revo- 
lution of the ISO satellite on 11 December 1995. The ISO 



satellite is described by Kessler et al. ( 1996 ). SWS and its ob - 
serving templates are described by de Graauw et al. ( 1996 ). 
The observations used the SWSOl template: a spectral scan 
from 2.4 to 45 /im. The observations were done at the slowest 
speed, which reduces the spectral resolution to roughly half the 
nominal resolution of SWS. The spectra were reduced with the 
SWS Interactive Analysis software (pipeline V4.3) and were 
subsequently smoothed with a gaussian profile of a half-width 
corresponding to a resolving power of 800. The wavelength 
calibration and the flux calibration were done with the set of 
calibration files adopted as a standard for this issue of A&A 
letters. The wavelength calibration is discussed by Valentijn et 
al. ( |1996| ); the flux calibration by Schaeidt et al. ( |1996[ ). The 
accuracy of the current flux calibration is estimated to be 30 %. 
Especially the SWS spectrum of NGC 7027 is extremely 
rich. More than 100 emission lines were identified, down to a 
flux of 2 X lO""" W/m^. It is impossible to discuss or even 
present all these features. This will be postponed to a later 
publication. In this paper we report the first detection of the 
[Arvi] 4.53 /im and the [Ne VI] 7.65 /im lines in the spectrum 
of NGC 7027. This makes Ne^+ the ion with the highest ioniza- 
tion potential observed in NGC 7027. In Table |l] the fluxes for 
these lines are compared with the 3 a upper Umits determined by 
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Table 1. Comparison of the observed line fluxes in the SWS spectra of 
NGC 7027 and NGC 6543 with previous observations in the literature. 
Where appropriate 3(7 upper limits are given. 



[Arvi] 



[Nevi] 



ion A 


SWS Lit. 
(lO-'-'W/m^) 


ref 


NGC 7027 



Bra 

[Mgiv] 
[ArVl] 
[Mgv] 
Pfa 

fNeVll 



4.05 
4.49 
4.53 
5.61 

7.46 
7.65 



9.3 
13.1 
10.6 
34.3 

24 



16.3 
31.0 

<2.0 

56.0 

5.6 

<5.4 



NGC 6543 



[Nam] 
[SIV] 
[Nelll] 
[OIV] 



7.32 
10.51 
15.55 
25.89_ 



c) Shureetal. (1983), 



0.23 
29.3 

71.1 
< 0.021 



43.0 

53.0 

<6.6 



a) Beckwith et afi 1984 ), b) Pottasch et al. ( 1986 ). 



Beckwith et al. ( |1984| ; hereafter B84) based on KAO spectra ob- 
tained in June 1981. The SWS spectra are shown in Fig. [I] We 
can see that the SWS detections yield a higher flux than the up- 
per limits given by B84. The first question to ask is whether flux 
could have been missed outside the aperture or if this discrep- 
ancy could be due to problems in the flux calibration in either of 
the spectra. The ionized part of NGC 7027 measures 9" x 13.5" 
based on the VLA 6 cm radio image of Masson ( 1989| ). This 
area fits completely in the aperture used by B84 (28" circular) 
and the smallest SWS aperture (14"x20"). However, in SWS 
observations pointing errors of up to 4" may occur The ab- 
sence of wavelength shifts indicates that any errors larger than 
about 2" could only have occuiTed along the slit (20") and we 
estimate that signal losses due to mispointing are less than 5 %. 

We also compared the fluxes of the four lines that were de- 
tected by B84 with the SWS fluxes (Table |I]). The average ratio 
of the B84 fluxes over the SWS fluxes is 2.0 ± 0.3. This might 
indicate a difference in flux calibration between the spectra. 
However, if anything, this would only exacerbate the discrep- 
ancy between the upper limits derived in 1 98 1 and the detections 
in 1996. We therefore feel confident that at least the [Ar VI] line 
and probably also the [NeVi] line have increased in strength 
over this period. We will discuss this further in Sect. Q 

We also report the non-detection of the [O IV] 25.9 fim line 
and the first detection of the [Nam] 7.32 /im line in the spectrum 
of NGC 6543. The upper limit for the [O IV] fine derived from 
the SWS spectrum is much lower than was previously measured 

|i. 



by Shure et al. (1983 ), as is shown in Table 



of this non-detection are discussed in Sect. 5 



The implications 



3. The photo-ionization models 

To model the planetary nebulae we used a modified version of 
the photo-ionization code CLOUDY 84.12a (Ferland, |l993| ). We 
assumed that the central star has a blackbody spectrum; that 
the nebula is spherically symmetric; that the density is constant 
inside the Stromgren sphere of the nebula, and varies as 1 /r^ 
outside; that the dust grains are intermixed with the gas at a 



100 


fl 


80 






- 


S 60 






: 










> 








'^ 40 




) 


20 


\ y 




y ^ 




, 


, ■ 



250 



200 




4.50 4.52 4.54 4.56 
^(|im) 



7.62 7.64 7.66 7.68 7.70 
^ (p.m) 



Fig. 1. The observed [ArVl] and [NeVi] lines in the spectrum of 
NGC 7027. 



constant dust-to-gas ratio and that the filling factor, describing 
the small scale dumpiness, is unity. 

For NGC 7027 the following observations were added to the 
SWS spectrum to constrain the modeling. The ultraviolet and 
optica l spectrum of NGC 7027 were taken from Middlemass 
( 199C ) together with the absolut e H/? fl ux. The distance to the 
nebula was taken from Pottasch ( 1996 ). The SWS spectrum of 
NGC 6543 was augmented with the ultraviolet spectrum from 
Middlemass et al. (1989) and the optical spectrum from AUer 
& Czyzak (1979). The absolute H/3 flux was taken from Acker 
et al. (1992) and the distance to the nebula from Van de Steene 



& Zijlstra (1994). Additionally also infrared continuum fluxes, 
radio fluxes and an angular diameter were used to constrain the 
models for both nebulae. 

To derive the physical parameters of the planetary nebulae 
from these observables, we adjust the model parameters until 
an optimal fit to the observables is found. For this we calculate 
a goodness-of-fit estimator x^, which is minimized. A detailed 
discussion of the method is given in van Hoof & Van de Steene 
(|1996|). 



Given the fact that in NGC 7027 the [Arvi] and Ukely also 
the [Ne Vl] line have increased in strength since 1981 one might 
expect that the optical and ultraviolet spectra have changed as 
well. To minimize the inconsistency we decided to omit the 
highest ionization levels from the older optical and UV spectra, 
because these are expected to be influenced most by an increase 
in stellar temperature. 

For NGC 6543 we made two models. The input for both 
models was identical except that in model 2b a modified black- 
body spectrum was used. This is discussed further in Sect. ^. It 
is well known that NGC 6543 has a very extended halo which 
contains most of the mass (Middlemass et al., 1989| ; Manchado 
&Pottasch, 1989). However, the assumption of constant density 
inside the Stromgren sphere does not permit us to model this 
halo. Since the line flux is expected to be completely dominated 
by the core region, this poses no problem. 

The results of the modeling are shown in Table g. Most 
symbols have a commonly accepted meaning. F denotes the 
dust-to-gas mass ratio and e the logarithmic abundance of an 
element (e(H) = 12.00). The quoted number density of hy- 
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Table 2. Parameters for the CLOUDY models of NGC 7027 and 
NGC 6543. Model 1 is for NGC 7027. Models 2a and 2b are for 
NGC 6543; these are identical except that different central star spec- 
tra were assumed. Abundances based on only one line are marked 
uncertain. 





model 1 


model 2a 


model 2b 


Teff (kK) 


161.4 


42.3 


58.1 


L, (Lq) 


6894. 


5023. 


3172. 


Tin (rape) 


2.4 


27. 


26. 


rstr (rape) 


22. 


55. 


52. 


rout (rape) 


345. 


55. 


52. 


log(nH) (cm-3) 


4.49 


3.56 


3.60 


Te(kK) 


16.13 


9.36 


7.57 


log(ne) (cra-3) 


4.55 


3.61 


3.64 


log(r) 


-2.52 


-2.20 


-2.06 


Mion (M0) 


0.058 


0.085 


0.070 


M,h (Mq) 


2.10 


0.085 


0.070 


e(He) 


10.98 


11.11 


10.94 


<C) 


8.87 


8.50 


8.96 


e(N) 


8.12 


7.65 


8.10 


e(0) 


8.65 


8.44 


9.09 


e(Ne) 


7.79 


7.81 


8.13 


e(Na) 


5.94 


6.31: 


5.82: 


e(Mg) 


6.82 






6(A1) 


4.87: 






e(Si) 


6.90: 


6.71: 


6.39: 


e(S) 


6.76 


6.91 


7.04 


e(Ar) 


6.34 


6.41 


6.49 


e(Ca) 


4.64: 






e(Fe) 


5.87 






DCpc) 


790. 


1080. 


1080. 


X^ 


5.4 


29.8 


8.6 



drogen is the value inside the Stromgren sphere. The electron 
temperature and density are averaged over the nebula. 



4. Discussion of the NGC 7027 results 

Since infrared fine-structure lines are rather insensitive to elec- 
tron temperature, the increase in strength of the [Ar VI] line must 
stem from a change in the ionization structure of the nebula. 
The ionization structure is the result of photo-ionization by the 
stellar radiation field and possibly also from shock heating due 
to the fast wind from the central star Such a wind has never been 
observed for NGC 7027, therefore the presence of very strong 
shocks can be ruled out. Our model predicts that the [Ar VI] line 
is formed in a region ranging roughly from 6 to 16 mpc from 
the central star. Based on the [Arvi] line strength, an Einstein 
A coefficient of A = 0.0966 cm~' (Mendoza, 1983) and assum- 
ing a relative upper level population of 1 % we can derive the 
minimum width of the line emitting region to be 1.7 mpc. As- 
suming 20 kK for the average electron temperature, the sonic 
travel time through the line emitting region would be > 70 yr, 
considerably longer than the 15 yr over which the spectrum 
changed. The electron temperature only enters as the square 
root in the sonic velocity, therefore it would have to be much 
higher than 20 kK to account for the discrepancy. This can be 
ruled out based on our models. Adding extra heating to the 
nebula (due to energy deposition by the stellar wind) equivalent 
to 10 % of all energy radiated by the central star would only 



double the electron temperature. Therefore we think that shocks 
are ruled out as the source of the variability in the spectrum. 

By changing the parameters in our photo-ionization model 
we could assess that changes in the nebular density, inner radius 
and dust-to-gas ratio have little effect on the line-strength of the 
[Arvi] line. An increase in the luminosity of the central star 
would increase the strength of all lines, which is not observed. 
Therefore the increased strength of the [Arvi] line must in 
all likelihood be attributed to a change in the spectral energy 
distribution of the central star, possibly an increase in effective 
temperature. The nebula reacts rapidly to a hardening of the 
stellar spectrum. We calculated the typical timescale needed to 
ionize Ar"^+ to be 0.2 yr. 

When we change the stellar temperature in our photo- 
ionization model, we see that the strength of the [Ar Vl] line 
indeed is sensitive to this parameter However, in order to ac- 
count for the large rise in strength of the [Arvi] line an un- 
realistic large increase in the stellar temperature of roughly 
60 kK to 100 kK is needed. Other lines show an increase in 
strength which is consistent with a more moderate increase in 
stellar temperature, or even with a decrease. These mixed re- 
sults might be caused by calibration problems since the changes 
for these lines are much more moderate. The evidence gathered 
so far does favor a change in the spectral energy distribution of 
the central star, possibly an increase in the stellar temperature. 
However, more observations are needed to confirm this result. 

5. Discussion of the NGC 6543 results 

In Table H we reported a stringent upper limit for the strength 
of the [Oiv] line in NGC 6543. This is not in contradiction 
with the detection of Oiv 1342 A P Cygni profile in the lUE 
spectrum of NGC 6543 since this UV line is formed in the wind 
coming from the central star, i.e. is formed much closer to the 
central star (Castor et al., 1981; hereafter CLS). We also report 
the detection of [Siv], [Nelll] and [Nam] emission. These 
detections, together with the upper limit, allow us to derive 
clear constraints on the shape of the central star spectrum. We 
will discuss this by starting off with ablackbody approximation, 
an assumption still often made in planetary nebula modeling. 

The ionization threshold is 34.8 eV for S^+, 41.0 eV for 
Ne2+, 47.3 eV for Na2+ and 54.9 eV for 0^+. The presence of 
the [Na III] line indicates that sufficient photons with energies 
above 47.3 eV are present to produce appreciable amounts of 
Na^+. When using a blackbody spectrum this makes it very 
likely that there are also enough photons to produce a significant 
amount of 0^+. The [O IV] line is expected to be a more sensitive 
tracer than the [Nam] line since the oxygen abundance is so 
much higher that, in all probability, it would more than make up 
for the lower transition probability of the [O IV] line. Therefore, 
if significant amounts of 0^+ were present in the nebula, the 
[O IV] line should be clearly detectable. From its absence we 
conclude that the blackbody approximation can not be valid; 
the stellar flux must make a considerable drop at the He II limit. 

This is also confirmed by our modeling. When we force the 
model to obey the upper limit on the [O IV] line when using a 
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Table 3. Comparison of the observed fluxes in NGC 6543 with model 
predictions. 



ion 


A 
(/im) 


SWS model 2a model 2b 




(10"'*W/m-) 


[Siv] 


10.51 


29.3 5.0 21.0 


[Nelll] 


15.55 


71.1 26.0 71.2 


[OIV] 


25.89 


< 0.021 0.027 0.0040 



blackbody spectrum (model 2a), we see that the excitation of 
the model spectrum is much lower than what is observed (see 
Table g). In general the fit is bad as is expressed by the high 
X^. The effective temperature is in good agreement with the 
measured H I Zanstra temperature of 47.0 kK (Kaler & Jacoby, 



1991). However, the predicted Hi Zanstra temperature based 



on the model spectrum is 31.0 kK, much lower than observed. 
We decided to test our hypothesis by modifying the black- 
body spectrum such that beyond the H I limit the flux would be 
multiplied by a constant factor 0.40 and beyond the He II limit 
by an additional factor 10^^. This resulted in model 2b. We can 
see that the overall fit of model 2b is much better and also that it 
has the right degree of excitation. The high oxygen abundance 
we find does not invalidate our model since the oxygen abun- 
dance in NGC 6543 is known to be high (Middlemass et al.. 



1989; Manchado & Pottasch 1989) 



The effective temperature we find is higher than the upper 
limit for the stellar temperature derived by CLS. However, we 
argue that this upper limit is not valid. It is based on the He II 
Zanstra temperature derived from an upper limit for the nebular 
He n A 4686 emission. The Zanstra method makes the assump- 
tion of a blackbody stellar radiation field, which is clearly not 
valid. A stellar radiation field of any given effective temper- 
ature can satisfy the constraint on the nebular He II emission 
provided that the flux in the Hell Lyman continuum is low 
enough. This is also confirmed by our model. The predicted 
nebular Hell A 1640 and A 4686 emission is 2.25 x 10"^'' 
W/m^ and 2.46 x 10^'^ W/m^ respectively, well within the 
restrictions imposed by CLS. The predicted Hi Zanstra tem- 
perature based on the model spectrum is 41 .6 kK, in agreement 
with the observations. The stellar temperature has also been de- 
rived by Lucy & Perinotto ( 1987 1. Their preferred value of 5 1 .6 
kK is higher than the value adopted by CLS, but still lower than 
the value we find. Our data point would fall to the left of the 
curve in their Fig. 2. The effective temperature is sensitive to the 
adopted multiplication factor at the H I limit. Since the nebular 
spectrum is not particularly sensitive to this factor, we decided 
to fix the value such that the predicted stellar V ma gnitude 
would coincide with the value from Castor et al. (1981 ). 



Photo-ionization modeling of the [O IV] line can prove to 
be a valuable test for predictions of stellar atmosphere models, 
since the results appear to be sensitive rather directly to the 
magnitude of the drop at the He II limit. This line is more suitable 
than the He II lines since there can be no ambiguity between 
stellar and nebular emission. 

Using the parameters of model 2b we determined that the 
stellar flux drops at least by a factor 350 at the He II limit in 



order to obey the upper limit on the [O IV] line. In a future paper 
we hope to present a more stringent upper limit. 

6. Conclusions 

In this letter we reported the first detection of the [Arvi] 
4.53 /im and [Nevi] 7.65 /xm lines in the spectrum of 
NGC 7027. The strength of the [Arvi] line and likely also 
the [Ne VI] line have increased since 1981. The most likely ex- 
planation for this variability is a change in the spectral energy 
distribution of the central star, possibly an increase in effec- 
tive temperature. However, further observations are needed to 
confirm this result. 

We also reported a non-detection of the [O IV] 25 .9 /im line, 
and the first detection of the [Nam] 7.32 /im line in the spectrum 
of NGC 6543. The ionization energy needed to produce 0^+ is 
just beyond the He II limit, and the absence of this line shows 
that the stellar flux drops at least by a factor 350 at the He II 
limit. Modeling the [O IV] line may prove to be a valuable test 
for atmosphere models. 
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